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Abstract. Motivated by the string landscape, the inflationary potential may be high dimen¬ 
sional and complicated. We propose a new way to construct some high dimensional random 
potentials, and simulate inflationary dynamics on top of that, for up to 50-dimensions in 
field space. Especially, random bifurcations of classical inflationary trajectory are studied. 
It is shown that the bifurcation probability increases as a function of number of dimensions. 
Those random bifurcations are not consistent with observations, and dramatically limit the 
parameter space of inflation on a complicated landscape. Eor example, in 10 dimensions, 
only 10“^ ~ 10“® of the parameter space volume leads to unique classical trajectories. The 
rest is ruled out by random bifurcations. 
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1 Introduction 

Inflation is the leading paradigm for the very early universe cosmology. However, the de¬ 
tailed dynamics of inflation still remains mysterious. We have not confirmed whether the 
inflationary potential is single or multiple dimensional, smooth or bumpy, featuring straight 
trajectory or curved, etc. 

The situation is rapidly improving with the on-going and future experiments. With 
more data on the cosmic microwave background (CMB) and the large scale structure (LSS), 
we expect to have better confidence in addressing the above questions. On the other hand, 
it is also important to better understand the underlying theoretical construction of inflation, 
so that the naturalness of those models can be better addressed. 

One of the particularly interesting ideas for theoretically constructing inflationary mod¬ 
els is the string landscape. The string landscape predicts that string theory has very compli¬ 
cated vacuum structure. A large number of light fields can be present in the string landscape, 
unless stabilized explicitly. The number of those fields can easily reach 0(100). Inflation on 
such a random potential has been studied in [1-8]. 

To study inflation on a random potential, the first step would be to construct such 
a random potential numerically There have been a few methods for random potential 
generation in the literature. One can either generate random numbers in position space 
and interpolate [9]; or generate Fourier coefficients and then Fourier transform to field space 
[10, 11]; or making use of Dyson Brownian motion [12, 14]; or use Gaussian covariance 
function to generate the potential dynamically [13]. Those methods, looking from a different 
angle, can also be classified into two types: 

• Static methods: A whole patch of inflationary potential is generated (or at least deter¬ 
mined before evolving the equation of motion of the inflaton. The pros and cons of 
static methods of potential generation include 

^There are analytical studies of statistical features of such random potentials. Even in those case, it would 
be helpful to numerically test the analytical calculation, and extend it where perturbative methods break 
down, with numerical methods. 

^It may be interesting to investigate the possibility to determine, instead of actually generate the whole 
potential. The time and space complexities may be greatly reduced. But such algorithms are not available, 
to the best of our knowledge. 
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— The potential is non-perturbative in the sense that we are not Taylor expanding 
the potential around any point. 

— There is no bias from different classical trajectories of inflation, because the po¬ 
tential is trajectory-independent. 

— Scenarios such as bifurcation, or self-crossing of inflationary trajectory can be 
studied. 

— Unfortunately, it is computationally extremely expansive to generate such a po¬ 
tential in high dimensional field space. Because the computational complexity to 
generate such a potential is 0{e^). To calculate such a fully random potential, 
the exponential computational expense is unavoidable because even to generate 
random numbers on the vertices of an N dimensional cube needs 2^ random num¬ 
bers. Interpolation of those numbers to obtain continuous fields is another 0{e^) 
calculation. 

• Dynamical methods: The potential is dynamically generated along an inflationary tra¬ 
jectory. The advantages and disadvantages of dynamical methods are to the opposite 
compared with static ones: 

— Fast computation for the existence of second derivative). 

— The potential only applies for a neighbourhood of the inflationary trajectory, be¬ 
cause the dynamical methods depend on Taylor expanding around the inflationary 
trajectory. 

— It is hard to study scenarios such as bifurcation, or self-crossing of inflationary 
trajectory. For example, if the trajectory is self-crossing, the potential at the self¬ 
crossing point shall have different values and derivatives when the inffaton runs 
into this point from different trajectories. On the other hand, it is interesting 
to note that after taking special care for the consistency of crossing points, it is 
possible to study such behaviors using dynamical methods [13]. 

In the present paper, we investigate bifurcations in a high-dimensional random potential. 
In section 2, we propose a new approximation method of potential generation, combining the 
advantages of static and dynamical methods. 

After the construction of the random potential, we move on to study bifurcations. When 
hitting bumps in the potential, the inflationary trajectory may bifurcate. The situation is 
illustrated in Fig. 1. The situation is known as multi-stream inflation [15]. 

Bifurcations of inflationary trajectory, if under control (for example, follows from spon¬ 
taneous symmetry breaking), can provide explanations of CMB statistical asymmetries [16] 
and/or the CMB cold spot [17]. This is because the e-folding number difference between 
different patches of the universe are slightly different. Moreover, temporary domain walls 
form between different trajectories before those trajectories recombine, which contributes to 
CMB statistical anisotropies [18]. The dynamics of eternal inflation can also be affected by 
such bifurcations [19, 20]. 

However, if the bifurcations happens during observable inflation, while not under con¬ 
trol, as in the case here on a random potential, disasters happen [9]. This is because random 
bifurcations typically have e-folding number difference of order 1 in different branches. Thus 
the CMB, observed at different directions, shall look highly different from each other, with 
AT/T ~ 0(1). This contradicts the CMB/LSS observations on large scales, and primordial 
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Figure 1. Multi-stream inflation with bifurcations of classical inflationary trajectory. 


black holes form after such scales return to the horizon, making contradicting predictions for 
small scale physics all the way down to reheating. 

Thus we should not allow random bifurcations to happen. As a result, inflation on 
a random potential is constrained. In [9], the bifurcation constraint on a two dimensional 
potential is considered. We find that the constraint is pretty weak for the case of a two 
dimensional potential, if the characteristic lengths of the random features on the inflation 
and isocurvature directions are equal (i.e. the random features are statistically circular 
shaped). However, if there is an ellipticity, such that the random feature is statistically more 
than 10 times longer in the inflation direction, bifurcations can happen easily. The situation 
is illustrated in Fig. 2. Examples of such potential with ellipticity can be found in [21]. 



Figure 2. Relation between ellipticity and bifurcation probability. 
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In this paper, we shall keep the features statistically circular, and consider higher di¬ 
mensional random potentials. The method of calculation is similar to that of [9], which is 
reviewed in Section 3. In Section 4, we present the result of bifurcation constraint up to 
50 field dimensions, for large and small field inflation models respectively. We conclude in 
Section 5. 

2 Generation of random potential 

In this section, we propose a method for high-dimensional random potential generation. Our 
purpose is to combine the advantages of static and dynamical methods, and enable study of 
bifurcation on a high dimensional random potential. 

We shall use 


V = Vo{(l>i) + Y,5yi,MiAj) , (2.1) 

i<j 

where (pi is the inflation direction with a slope, and the other directions are flat up to random 
noises 6Vij{(f>i, <pj). 

If we were to use such a multiple dimensional potential for a general purpose, there 
shall be serve draw-backs, because the potential does not have enough “randomness”. For 
example, consider 3-dimensional case. We have 

F = Fo(<^l) + 6Vi,2{(pl,cP2) + SVi,3{(Pi, 03) + <51/2,3(02, 03) • (2.2) 

We observe that = 0, which is not true for a fully random potential. Also, 

depends on 0i, which is also unpleasant. 

However, for the case of inflation with one dominant direction, such draw-backs shall 
not cause problems for practical purposes. Because for the purpose of studying bifurcations, 
we shall only be interested in the first derivative of the potential. (And the potential can 
easily be generated into a summation of l/i,j,fc(0i, 4>ji 0fc) or higher orders if needed.) 

Although the current method is a simplified version of static potential generation 
method, it is also worth mentioning the connection with the dynamical potential generate 
methods. During inflation, the inflaton field rolls much faster compared with other massless 
or stabilized directions This is because, at each e-fold, the inflaton rolls a distance of (p/H, 
plus small quantum fluctuations. On the other hand, the other directions moves a distance 
of because of quantum fluctuations. The ratio of those amplitudes is 1 /~ 10^. 

As a result, effectively, 0i acts as the time variable along the inflationary trajectory. 
Thus the method can be considered effectively as dynamical potential generation, where 0i 
parametrizes the inflationary trajectory. However, thanks to the static nature, the current 
method shall not suffer from the limitations of local Taylor expansion, and can be used to 
study bifurcation events. 

The time complexity for this potential generation method is 0{N‘^) if only random first 
derivative is needed (as in our case in later sections). It is worth to note that, the algorithm 
may further be optimized into 0{N), by noting the above fact that 0i parametrizes the 
inflationary trajectory. One can construct V = Vb(0i) -|- <5^i(</’i,'/’i)- However, if this 

^If there are a few classically rolling directions, one can always rotate the fields such that only one field rolls 
classically. Moreover, a turning of trajectory can be parametrized as interactions between different directions 
with straight trajectory from non-trivial field space metric [22, 23]. 
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simplification is made, statistical symmetry of field rotation in 5V is broken. As a result, 
if one still hope the random part of the potential to be symmetric, additional tunings are 
needed. To avoid possible issues, we shall not use this simplification in this paper. 

For doing numerical work, we yet need to specify a method to generate 6Vij{4>i, 4>j). As 
mentioned in the introduction, there are many ways to do it. Here we shall first generate 
on a grid, then interpolate to get smooth helds. 

3 Dynamics of inflation 

In this section, we review the inflationary dynamics. On a multi-dimensional landscape, the 
classical equation of motion for the fields is given by 


+ diV — 0 , (3-1) 

where the potential V is constructed in the previous section. The shorthand notation di is 
used to denote the partial derivative with respect to the held (pi. To simplify the calculation, 
H ~ (SM^) is still a good approximation. This is because, on a random potential, in 

order for the inhationary trajectory not to get stuck, the kinetic energy of the inhaton is not 
likely to exceed twice the slow roll value, thus still a few percents. 

To encode quantum huctuations, we use Starobinsky’s stochastic method [24], where 
quantum huctuations are considered effectively as source terms for the classical equation of 
motion. Such source terms can be derived by integrating out the sub-Hubble modes. The 
quantum-corrected equation of motion reads 

Q 

Pi + 3HPi + d^V = —Hirjiit) , (3.2) 

ZTT 

where the rji term is the random source to denote the quantum corrections and follows 
independent Gaussian distribution The Tji terms are normalized as 

{Vi{t)r]j{t')) = 5{t - t')6ij , (3.3) 

such that during a Hubble time and averaged on a Hubble volume, the amplitude of quantum 
huctuation on each held direction is H/{2'k). 

In the numerical calculation, we have to discretize the time interval as At = tn — tn-i- 
While replacing differential equations into difference equations, the Dirac delta function is 
replaced by 


^ita — tb) —)• -—— , (3.4) 

where tn denotes the time of the nth step. 

In the following section, inhation on a random potential is studied based on the above 
setups. 

"^For the isocurvature fields, the pi could be small, or even zero without the presence of the RHS. In such 
a special case, the RHS is the major driving force for this differential equation. (But note that this can never 
happen in the infiaton direction.) However, this fact does not indicate that the back-reactions are out of 
control. This is because, the scalar field fluctuation amplitude H/(2k) originates from de Sitter expansion. 
Thus as long as H does not receive large back-reaction (which is always true in our work), the back-reaction 
to Eq. (3.2) is small and the random source is highly Gaussian. 
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4 Bifurcations on a random potential 

In this section, we investigate the bifurcation probabilities with different field-space dimen¬ 
sions. As we shall show, the bifurcation probabilities increase as a function of field-space 
dimension, in both large field and small held inhation models. 

4.1 Large field inflation 

First, we consider large held inhation, with potential Vo((/>i) = , plus a random source 

introduced following the algorithm described in Sec. 2. The mass parameter is chosen to be 
m = 6.5 X Note that we have picked the ijii direction to be the direction with a 

non-random slope. 

Parameters are introduced to parametrize the random part of the potential: 

• We use Ap(/)j (i = l,...,n) to denote the characteristic distance in the held space, 
between random bumps in the potential. Operationally, we are generating the ran¬ 
dom potential by interpolation on grids. Thus Ap(/>j becomes the grid size in the f-th 
direction. 

In the two dimensional example [9] , it is shown that the bifurcation probability increases 
signihcantly as a function of Ap4>i/Ap(f)2. This is because, when Ap4>i increases, bumps 
in the cpi direction are more distant from each other. Thus it is easier to keep the 
inflationary trajectory slowly rolling. On the other hand, when Ap(p 2 decreases, bumps 
in the (j )2 direction are narrower. Thus it is easier for a small quantum fluctuation in 
the (/)2 direction to bifurcate the inflationary trajectory. We expect that this intuition 
still hold in general n-dimensions, simply by replacing Ap(j )2 into ApCpj (j = 2, ...n). 

In the numerical study, to reduce the number of parameters, we shall choose Ap4>i = 

Ap(f>2 — ... — ApCpn — A. 

• We use Ai (z = I,..., n) to denote the relative amplitude between the random part of 
the potential and the slope Vq, contributed from the f-th dimension. For simplicity, we 
calculate only the case Ai = ... = An = A. As A increases, the parameter space are 
parted into three regions: 

— Small A\ If there is no randomness {A = 0), we essentially have single held 
inflation, plus decoupled isocurvature fluctuations. Thus it is natural to expect 
that for sufficiently small A, no bifurcation or stuck behavior takes place. 

— Intermediate A: As we shall show, a bifurcation region emerges as A increases. 
This region is largely overlooked in the previous study of multi-field inflation 
(except [9]), and should be excluded, because bifurcation introduces too large 
density fluctuations. 

— Large A: When A further increases, the slope diV shall not be monotone. Thus 
the inflationary trajectory gets easily stuck at eternal inflation. We shall exclude 
this case from our study. 

It is intuitive to think about bifurcations from Fig. 1. Nevertheless, it is important to 
define the bifurcation events more precisely. We define bifurcation as that, at least one of 
direction cjij [j = 1, ...n) satisfies the condition \(j)j — > A at the end of inflation, where 

and (f)^ are two sample trajectories. Two further comments are in order for the definition 
of bifurcation probability 
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• We have excluded the eternal inflation samples when defining the bifurcation probabil¬ 
ity. 

• The above definition of bifurcation stops making sense when A < y/Wg,H, where W ~ 
60 is the e-folding number of observable inflation. This is because when A < 
trivial random fluctuations may be identified as an bifurcation event following the above 
definition. However, the parameter space that we have considered, in both large field 
and small field examples, A is orders of magnitude greater than \/7^H. Thus this 
subtlety does not affect our study. 

In the Fig. 3 we study the bifurcation probability numerically for large field inflation. 
The pink region is the constraint by stuck or unphysical ending of our numerical calculation 
(in which case the field runs too far away in the isocurvature direction, beyond the size of 
the grid) , which is beyond the scope of our interest. 

Note that in the 3, 5 and 10 dimensional cases in Fig. 3, the stuck region (pink) does 
not fully cover the parameter space from above. To understand this, it is helpful to take a 
closer look at how the trajectory can get stuck in a highly random potential. To get stuck, 
the inflaton typically hrst falls into a potential well from one side, and then cannot climb 
out from the other side. Note that the inflaton first obtains some kinetic energy when it falls 
into a potential well. After that, there are two cases when the inflaton get stuck and cannot 
get out: 

1. The potential well at “the other side” is higher than the potential well at “one side” 
plus the inflaton kinetic energy (calculated before falling into the potential well). In 
this case, there is no way for the inflaton to classically climb up (and the quantum 
effect is small) from energy conservation. 

2. The potential well at “the other side” is higher than the inflaton kinetic energy; but 
not as large as the kinetic energy (calculated before falling into the potential well) plus 
the energy obtained from falling into the potential well from the “one side”. In this 
case, whether or not the inflaton gets stuck depends on the feature distance A. This is 
because, when A is small, the fall-climb happens very fast and the fraction force from 
Hubble parameter does not play a significant role. In this case, the inflaton can climb 
out following energy conservation. However, when A is large enough, the additional 
inflaton kinetic energy which is obtained from falling into the potential well gets diluted 
from Hubble expansion. As a result, the inflaton cannot climb out if the potential well 
is higher than the inflaton kinetic energy. 

Note that case 2 is less demanding to be triggered. Thus the stuck probability shrinks 
when the feature distance A is small. 

To better observe the scaling behavior as a function of number of fields, we fix A at 
O.lMp and consider a few examples of amplitudes. The result is shown in Fig. 4. One can 
clearly observe that the increment of bifurcation probability is faster than linear as field- 
space dimension increases (except that when the bifurcation probability is close to 1, the 
rapid increasing behavior certainly stops). This behavior is natural to expect. Because 
linear increment simply reflects the fact that there are more directions to have bifurcation 
events. On the other hand, bifurcation may happen in a (time dependent) linear combination 
of different directions, which further increase the bifurcation probability. As one can read 
from the figure, the increment behavior is approximately quadratic. 
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Figure 3. Density plot for bifurcation probabilities in 2, 3, 5 and 10 dimensions. The pink regions 
corresponds to exclusion regions where the inflationary trajectory is stuck at eternal inflation. Notice 
that in the 10-dimensional case, as the stuck region gets bigger, we have zoomed in the y-axis. Here 
the unit of A is Planck mass Mp. 


Finally, let us estimate the fraction of non-eternally-inflating parameter space which 
bifurcates. To actually calculate such a fraction of parameter space requires exponentially 
computational complexity. Because we need to All the n-dimensional parameter space Ai 
(i = 0,1,..., n) with numbers, where Ai can now take different values for different i. However, 
one can nevertheless estimate the parameter space by noting that the non-bifurcating region is 
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Figure 4. The bifurcation probability as a function of number of field-space dimensions, for large 
field inflation. 


a ball-like sub-volume of the non-eternally-inflating parameter space. Thus the volume ratio 
is exponentially small. The situation is illustrated in Fig. 5. In general, for N dimensions, 
the volume fraction / scales as 


/ = 



(4.1) 


where is the amplitude boundary for bifurcation (where more than 10% of the trajectories 
bifurcate), and is the amplitude boundary for eternal inflation. 

Note that with multiple dimensions, eternal inflation become unlikely because with more 
directions to go, the field is less likely to get stuck. Thus the non-eternal parameter space 
tends to be larger than the exponential estimate. However, with multiple dimensions, the 
probability of bifurcation is increasing because every direction, and their linear combinations, 
may bifurcate. Thus the non-bifurcating parameter space tends to be smaller than the 
exponential estimate. As a result, it is safe to estimate that the non-bifurcating parameter 
space is an exponentially small subset of the non-eternally inflating parameter space, as a 
function of N. 

To put in real data, the estimate of bifurcation probability is plotted in Fig. 6. For 
example, in 10 dimensions, only about 10“® of the non-eternal inflation parameter space 
does not suffer from bifurcations. 
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eternal inflation boundary 




Figure 5. In high field-space dimensions, the parameter space of non-bifurcating region takes an 
exponentially small volume. Here the pink boundary corresponds to the parameter space of non¬ 
eternal inflation. The gray ball-like shape corresponds to the non-bifurcation region. Note that this 
figure is for illustration of the scaling rule, which is not calculated from real data. 



Figure 6. An estimate of the non-bifurcating to non-eternal inflating volume ratio for different 
dimensions for large field inflation. 

4.2 Small field inflation 

For small field inflation, we consider a brane inflation [25] potential 

H = yo(l - 4) , (4.2) 

n 

plus a small random potential. The parameters are taken according to [26], with /x = O.OlMp. 

The bifurcation probability as a function of number of field-space dimensions is plotted 
in Fig. 7, where A = 0.0005Mp is chosen. Similar to the large field case, a faster-than-linear 


- 10 - 












growth is spotted. On the other hand, bifurcation requires much smaller A in small field 
inflation than large field inflation, because the field is rolling more slowly, and thus more 
sensitive to features in the potential. The ratio of non-bifurcating to non-eternal inflationary 
parameter spaces is estimated in Fig. 8. It is noted that in small field inflation, the reduce of 
parameter space is not that dramatic as the case of large field inflation. Nevertheless, only 
about 10“^ of the total non-eternal inflationary parameter space is free from bifurcations. 



Figure 7. The bifurcation probability as a function of number of field-space dimensions, for small 
field inflation. 
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Figure 8. An estimate of the non-bifurcating to non-eternal inflating volume ratio for different 
dimensions for small field inflation. 

5 Conclusion and outlook 

To conclude, after proposing a new way to approximately construct a random multi-field 
inflationary potential, we extend our previous study of bifurcation phenomena into the case 
of higher field space dimensions. The situation turns out to be not a qualitative correction, 
but instead a remarkable quantitative change. There are two sources of enhancements of the 
bifurcation probability: 

• Enhancement of bifurcation probability as a function of field-space dimensions, when 
fixing other potential parameters. As we have shown, when we add field directions, the 
bifurcation probability increases. The speed of increment is faster than linear in both 
the large field and the small field examples that we have studied. This result is natural 
to expect. This is because a linear behavior should be the lower bound of bifurcation 
probability as a function of field dimension, simply considering that there are more 
directions to move into. On the other hand, bifurcation could either happen on a linear 
combination of those field dimensions, which makes the behavior faster than linear. 

• A ball-like sub-volume in high dimensions has exponentially tiny space fraction com¬ 
pared to the whole parameter space. We do not know the exact volume ratio between 
the non-bifurcating sub-volume and the full parameter space. Because one has to sim¬ 
ulate a exponentially high number of potentials with different parameters to really 
account for this ratio. However, the ratio should scale exponentially as a number of di¬ 
mensions. Thus only an exponentially small fraction of parameter space is left without 
bifurcations, which has a chance to agree with any observations. 

In the present paper, we have not applied constraints from observation, like features in 
the primordial power spectrum, and non-Gaussianities, to limit the shape of the inflationary 
potential. Other than technical simplicity, we have the following considerations for not mix-in 
the latest observations: 
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• The current observational constraints on highly scale dependent features are not yet 
very strong. Especially, the current data-fitting typically focuses on the encounter of a 
single feature instead of averaging a large number of features. 

• Most of the experiments, especially the most precise ones, observes only the first 10 e- 
folds of inflation, among the 50-60 e-folds along the “observable” inflationary trajectory. 
On the other hand, our bifurcation constraint covers a much greater number of e-folds 
because otherwise primordial black holes could have been formed from the density 
fluctuation of the bifurcated trajectory (but those primordial black holes cannot be 
predicted if one followed only one classical trajectory), which contradicts observations. 

• In numerical studies of observational features of random potentials (such as non- 
Gaussianities), one typically first generate a random potential, and test that the po¬ 
tential does not lead to eternal inflation. Then a statistical study of the observational 
features is performed (see, for example, [5]). Our result show that, an exponentially 
large portion of parameter space should be cut off well before hitting the eternal infla¬ 
tion constraint. 

On the other hand, it is indeed interesting to combine all observational constraints and we 
hope to do so in future work. 
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